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ABSTRACT 

We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright 
continuum sources. We have used the GREAT instrument on SOFIA to observe the SH 1383 GHz ^n 3/2 7 = 5/2 «— 3/2 lambda 
doublet towards the star-forming regions W31C, G29.96-0.02, G34.3-I-0.1, W49N and W51, detecting foreground absorption towards 
all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H 2 S lio- loi (169 GHz), CS 7 = 2 - 1 
(98 GHz) and SO 82 - 2i (99 GHz) transitions. Upper limits on the H 3 S^ lo - Oo (293 GHz) transitions were also obtained at the 
IRAM 30 m. In nine foreground absorption components detected towards these sources, the inferred column densities of the four 
detected molecules showed relatively constant ratios, with A(SH)/A(H 2 S) in the range 1.1 - 3.0, A(CS)/A(H 2 S) in the range 0.32 - 
0.61, and A(SO)/A(H 2 S) in the range 0.08 - 0.30. The column densities of the sulphur-bearing molecules are very well correlated 
amongst themselves, moderately well correlated with CH (a surrogate tracer for H 2 ), and poorly correlated with atomic hydrogen. 
The observed SH/H 2 ratios - in the range 5 to 26 xl0“® - indicate that SH (and other sulphur-bearing molecules) account for 1% 
of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by 
standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by 
elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent 
dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM 
is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the 
column densities of H 2 S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity 
drift, relative to H 2 , of SH molecules produced by the dissociative recombination of H 3 S^ - reduce this discrepancy to a factor ~ 3. 
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1. Introduction 


To date, more than a dozen molecules containing sulphur have 
been detected unequivocally in interstellar gas clouds and/or 
circumstellar outflows: SH, SH^, CS, NS, SO, SiS, H 2 S, SO 2 , 
HCS+, OCS, C 2 S, C 3 S, HNCS, HSCN, H 2 CS, and CH 3 SH. All 
but one of these molecules were first detected using ground- 
based telescopes operating at millimeter or submillimeter wave¬ 
lengths, the only exception being the mercapto radical, SH, 
which has its lowest rotational transition at 1383 GHz, a fre¬ 
quency inaccessible from ground-based observatories and also 
lying in a gap between Bands 5 and 6 of the Herschel Space 
Observatory's HIFI instrument. The first detection of interstel¬ 
lar SH was obtained (Neufeld et al. 2012) in absorption towards 
the star-forming region W49N using the German Receiver for 
Astronomy at Terahertz Frequencies (GREAT; Heyminck et al. 


2012) instrument on the Stratospheric Observatory for Infrared 
Astronomy (SOFIA; Young et al. 2012). 

The chemistry of interstellar molecules containing sulphur is 
distinctive in that none of the species S, SH, S^, SH^, or H 2 S^ 
can react exothermically with H 2 in a hydrogen atom abstraction 
reaction (i.e. a reaction of the form X - 1 - H 2 —> XH - 1 - H.) This 
feature of the thermochemistry greatly inhibits the formation of 
sulphur-bearing molecules in gas at low temperature. 

In this paper, we present a corrected analysis of the SH abun¬ 
dance presented in Neufeld et al. (2012), together with addi¬ 
tional detections of SH obtained with SOFIA/GREAT toward 
four other bright continuum sources - W51, W31C, G34.3-1-0.1 
and G29.96-0.02 - and ancillary observations of H 2 S, CS, SO 
and H 3 S^ obtained with the IRAM 30 m telescope at Pico Veleta 
toward all five sources. (The observations of H 3 S^ led only to 
upper limits, not detections.) The observations and data reduc¬ 
tion are described in §2 below, and results presented in §3. The 

* GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA / Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DER Institut fiir Planetenforschung. 


1 




Neufeld et al.: Sulphur-bearing molecules 


Table 1. List of observed sources 


Source 

R.A. 

J2000 

Dec 

J2000 

Distance" 

kpc 

^sys 

km/s 

Date 

(SOFIA) 

SH 

Integration time (s) 
HjS CS SO 

H 3 S+ 

W31C 

18:10:28.70 

-19:55:50.0 

4.95 

-4 

17 Jul 2013'’ 

677 

8886 

1831 

3357 

3697 

G29.96-0.02 

18:46:03.80 

-02:39:20.0 

5.26 

98 

17 Jul 2013'’ 

752 

9213 

3682 

3682 

5531 

G34.3-t0.1 

18:53:18.70 

-h01:14:58.0 

3.8 

59 

22 Jul 2013'’ 

677 

9206 

5495 

5495 

3711 

W49N 

19:10:13.20 

+09:06:12.0 

11.11 

8 

28 Sep 201L 

495 

6400 

1221 

1221 

2737 

W51 

19:23:43.90 

+ 14:30:30.5 

5.1 

55 

02 Nov 2013" 

752 

6410 

3371 

3371 

3039 


References for distance estimates in order listed: Sanna et al. (2014); Zhang et al. (2014); Fish et al. (2003); Zhang et al. (2013); Sato et al. 
(2010). All are trigonometric estimates except that for G34.3-1-0.1, which is a kinematic estimate. 

® SOFIA deployment from Christchurch, New Zealand 
SOFIA deployment from Palmdale, CA, USA 


Table 2. List of observed spectral lines 


Molecule 

Transition 





Rest frequency 
GHz 

Beam size 

arcsec 

Ev/k 

K 

El Ik 

K 

^UL 

S“^ 

SH 

^Hj/z 7 = 5/2^ 3/2 

F = 

2+ 


2- 

1382.9056 

21 

66.4 

0.0 

0.47 X 10“^ 


2n3/2 J = 5I2^ 3/2 

F = 

3+ 


2- 

1382.9101 




4.72 X lO-'* 


2n3/2 7 = 5/2^ 3/2 

F = 

2+ 


L 

1382.9168 




4.24 X 10-3 


^Hj/z 7 = 5/2^ 3/2 

F = 

2- 

«- 

2+ 

1383.2365 




0.47 X 10-3 


2n3/2 7 = 5/2^3/2 

F = 

3 - 


2+ 

1383.2412 




4.72 X 10-3 


2n3/2 7 = 5/2^3/2 

F = 

2- 

<- 

1+ 

1383.2478 




4.24 X 10-3 

H2S 

Jk„.Kc = Lo ^ loi 





168.7628 

14.6 

27.88 

19.88 

2.68 X 10-3 

CS 

7 = 2 <- 1 





97.9810 

25.1 

7.05 

2.35 

1.67 X 10-3 

SO 

7a( = 32 2; 





99.2999 

25.6 

9.23 

4.46 

1.13 X 10-3 

H3S+ 

II 

p 

II 

T 

0 





293.4572 

8.4 

14.08 

0.00 

2.94 X 10-" 


abundances of SH, H 2 S, CS, and SO - and upper limits on the 
H 3 S^ abundances - derived in diffuse molecular clouds along 
the sight-lines to these bright continuum sources are discussed 
in §4 with reference to current astrochemical models. 

2. Observations and data reduction 

Table 1 lists the continuum sources toward which SH, H 2 S, CS, 
SO and HaS^ were observed, along with the exact positions tar¬ 
geted, the estimated distances to each source, the source sys¬ 
temic velocities relative to the Local Standard of Rest (LSR), the 
dates of each SOFIA observation, and the total integration time 
for each transition. In Table 2, we list the transitions targeted 
for each species, together with their rest frequencies, telescope 
beam size at the transition frequency, upper and lower state en¬ 
ergies, and spontaneous radiative decay rates. All the molecular 
data presented here were taken directly from the CDMS (Muller 
et al. 2005) or JPL (Pickett et al. 1998) line catalogs, with the ex¬ 
ception of those for the H 3 S^ transition; for the latter transition, 
which is not included in either catalog, we adopted the frequency 
measured by Tinti et al. (2006), and obtained the radiative decay 
rate for an assumed dipole moment of 1.731 D (Botschwina et 
al. 1986). 

All the SOFIA data reported here were obtained with the 
LI receiver of GREAT and the AFFTS backend; the latter pro¬ 
vides 8192 spectral channels with a spacing of 183.1 kHz. As 
described in Neufeld et al. (2012), the W49N SH data were ac¬ 
quired in the upper sideband of the LI receiver; observations 
of the other four sources were performed with the SH line fre¬ 
quency in the lower sideband. At this frequency, the telescope 


beam has a diameter of ~ 21" HPBW. The observations were 
performed in dual beam switch mode, with a chopper frequency 
of 1 Hz and the reference positions located 60" (75" for W49N) 
on either side of the source along an east-west axis. 

The raw data were calibrated to the (“forward beam 
brightness temperature”) scale, using an independent ht to the 
dry and the wet content of the atmospheric emission. Here, the 
assumed forward efficiency was 0.97 and the assumed main 
beam efficiency for the LI band was 0.67 (except in the case 
of the W49N data from 2011, where the beam efficiency was 
0.54). The uncertainty in the flux calibration is estimated to be 
~ 20% (Heyminck et al. 2012). Additional data reduction was 
performed using CLASlfl This entailed the removal of a base¬ 
line and the rebinning of the data into 10 -channel bins, resulting 
in a 0.40 km/s channel spacing for the final data product. 

All the IRAM data reported here were obtained using the 
EMIR receivers (Carter et al. 2012) on the 30-m telescope at Pico 
Veleta on 2014 March 16-18, with the H 2 S, CS, SO, and H 3 S+ 
transitions observed respectively using the El, EO, EO, and E3 
receivers. Since two receivers can be used simultaneously (and 
upper and lower sideband spectra obtained separately and simul¬ 
taneously for two orthogonal polarizations), we targeted the H 2 S 
transition in all observations together with either ( 1 ) the H 3 S^ 
transition when the weather conditions were excellent (precip- 
itable water vapor < 1 mm) or (2) the CS and SO transitions. The 
FTS backend was deployed for all observations, providing a 200 
kHz channel spacing over a 4.05 GHz bandwidth. In addition, 
the VESPA autocorrelator was used to target the CS transition 


^ Continuum and Line Analysis Single-dish Software. 
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Table 3. Continuum antenna temperatures and r.m.s. noise 


Source 

Line 


rA(cont) 

(K) 

T'rms 

(mK) 

Channel 

width 

(km/s) 

S/N'' 

ratio 

W31C 

SH 

1383.2 GHz 

4.11 

58.7 

0.40 

70 


H2S 

168.8 GHz 

0.75 

3.8 

0.35 

198 


CS 

98.0 GHz 

0.81 

10.8 

0.24 

75 


SO 

99.3 GHz 

0.78 

9.4 

0.24 

83 


H3S+ 

293.5 GHz 

0.86 

12.3 

0.20 

69 

G29.96-0.02 

SH 

1383.2 GHz 

1.47 

67.3 

0.40 

21 


H2S 

168.8 GHz 

0.36 

3.9 

0.35 

92 


CS 

98.0 GHz 

0.42 

6.4 

0.24 

65 


SO 

99.3 GHz 

0.43 

5.5 

0.59 

78 


H3S+ 

293.5 GHz 

0.41 

10.7 

0.20 

37 

G34.3+0.1 

SH 

1383.2 GHz 

4.38 

66.0 

0.40 

66 


H2S 

168.8 GHz 

1.47 

4.6 

0.35 

320 


CS 

98.0 GHz 

1.32 

5.9 

0.24 

222 


SO 

99.3 GHz 

1.35 

5.0 

0.59 

269 


H3S+ 

293.5 GHz 

1.53 

14.7 

0.20 

104 

W49N 

SH 

1383.2 GHz 

3.54 

208.3 

0.40 

17 


H2S 

168.8 GHz 

1.46 

5.3 

0.35 

278 


CS 

98.0 GHz 

1.57 

11.2 

0.24 

140 


SO 

99.3 GHz 

1.54 

12.7 

0.24 

120 


H3S+ 

293.5 GHz 

1.63 

16.9 

0.20 

96 

W5I 

SH 

1383.2 GHz 

4.82 

67.6 

0.40 

71 


H2S 

168.8 GHz 

0.96 

7.9 

0.35 

120 


CS 

98.0 GHz 

0.62 

6.9 

0.60 

90 


SO 

99.3 GHz 

0.62 

6.9 

0.59 

90 


H3S+ 

293.5 GHz 

1.41 

28.4 

0.20 

49 


Signal-to-noise ratio on the continuum 


toward W49N, W31C, G34.3-H0.1 and G29.96-0.02 and the SO 
transition toward W49N and W31C; the VESPA backend was 
used to obtain a 80 kHz channel spacing over a bandwidth of 
160 MHz. All the data were acquired in double beam switching 
mode, using the wobbling secondary with a throw of +45 arc- 
sec. Observations of Mars and Mercury were used for periodic 
checks of the telescope focus. 

The CLASS software package was used to coadd all the data 
obtained for each source, with a weighting inversely propor¬ 
tional to the square of the radiometric noise in each scan. The 
final coadded spectra show the antenna temperature, averaged 
over both polarizations, as a function of the frequency in the 
Local Standard of Rest (LSR) frame. 


3. Results 

In Table 3, we report the continuum brightness temperatures 
measured at each frequency, together with the r.m.s noise in the 
coadded spectra. The continuum spectra of H 2 S, SO and CS are 
of a particularly high signal-to-noise ratio, providing exquisite 
sensitivity to weak spectral features. In Figures 1 through 5, we 
show the spectra obtained toward each source. The left panels 
are scaled to show the full range of observed antenna temper¬ 
atures. The H 2 S, CS, and SO spectra show strong emission at 
the systemic velocity of each source, together with multiple nar¬ 
rower absorption features. The SH spectra show only absorption. 
The HaS^ spectra, which show multiple interloper emission lines 
(i.e. emission lines other than the targeted line), show no evi¬ 
dence for H 3 S^ emission or absorption. The right panels show 


the same data with the vertical scaling adjusted to reveal weaker 
absorption features. 

In each panel in Figures 1-5, the red parts of each histogram 
denote spectral regions that were assumed to be devoid of fore¬ 
ground absorption. These regions were subsequently ht with a 
flat continuum and one or more Gaussian emission features, to 
provide a background spectrum of the source itself, TbgCvLSR)- 
The observed spectra were then divided by that fitted back¬ 
ground spectrum (including emission lines), to obtain a trans¬ 
mission spectrum, f = T^IT^g, and an optical depth spectrum, 
t(vlsr) = -In f(vLSR)- In the case of the SH transition, which 
is split by lambda doubling and hyperfine splitting, the velocity 
scale is computed for the ^n 3/2 J - 5/2 «— 3/2 F - 3^ 2^ 

transition^ As is clear from Figures 1-5, the absorption pat¬ 
tern is therefore repeated with a velocity shift of +71.8 km/s, 
corresponding to the magnitude of the lambda doubling, and is 
broadened by the hyperfine structure. 

Figures A. 1 through A.5 (black histogram), in Appendix A, 
show the optical depth as a function of LSR velocity for the 
SH, H 2 S, CS, and SO transitions. The vertical axis shows the 
total molecular column density per velocity interval. The latter 
is computed for an excitation temperature of 2.73 K, the tem¬ 
perature of the cosmic microwave background. The validity of 


^ With the rest frequencies listed in Table 2, the velocity centroids of 
the SH absorption features were found to be systematically offset from 
those of other molecules. To bring them into accord, we adopted a +2 
MHz velocity shift for all the SH frequencies in Table 2; a correction of 
this magnitude is consistent with the experimental uncertainties in the 
laboratory spectroscopy. Figures 1-5 include this frequency correction. 
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W49N SH W49N SH 




W49N HjS 



LSR velocity (km/s) 
W49N CS 


W49N HjS 



W49N CS 




W49N SO 



W49N H,S* 



W49N SO 



W49N H,S* 



Fig. 1. Spectra obtained toward W49N. Red sections of each plot indicate spectral regions assumed devoid of absorption for the 
purpose of fitting a continuum with emission lines. The right panels are zoomed versions of the left panels, with the scaling chosen 
to reveal weak absorption features. Magenta lines indicate the lambda doubling and hyperfine splittings for SH. 
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W31C SH 



W31C H 2 S 





W31C HjS* 



LSR velocity (km/s) 


W31C SH 



W31C H 2 S 



W31C CS 



W31C SO 



W31C H 3 S* 



LSR velocity (km/s) 


Fig. 2. Spectra obtained toward W31C. Red sections of each plot indicate spectral regions assumed devoid of absorption for the 
purpose of fitting a continuum with emission lines. The right panels are zoomed versions of the left panels, with the scaling chosen 
to reveal weak absorption features. Magenta lines indicate the lambda doubling and hyperfine splittings for SH. 
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W51 SH W51 SH 




W51 H 2 S 



W51 HjS 



W51 CS 



W51 CS 



W51 SO 




W51 H 3 S* 



LSR velocity (km/s) 


W51 H 3 S* 



Fig. 3. Spectra obtained toward W51. Red sections of each plot indicate spectral regions assumed devoid of absorption for the 
purpose of fitting a continuum with emission lines. The right panels are zoomed versions of the left panels, with the scaling chosen 
to reveal weak absorption features. 
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G34.3H-0.1 SH G34.3H-0.1 SH 




8 

^ 6 

G34.3H-0.1 

H 2 S 

- 

: 

m" 4 



(/) 

2 



0 




-20 0 20 40 60 80 100 

LSR velocity (km/s) 


G34.3H-0.1 HjS 



G34.3H-0.1 CS 



G34.3H-0.1 CS 



G34.3+0.1 SO 




G34.3H-0.1 H^S* 



LSR velocity (km/s) 


G34.3+0.1 HjS* 



Fig. 4. Spectra obtained toward G34.3-1-0.1. Red sections of each plot indicate spectral regions assumed devoid of absorption for the 
purpose of fitting a continuum with emission lines. The right panels are zoomed versions of the left panels, with the scaling chosen 
to reveal weak absorption features. 
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2.0 

1.5 

^ 1-0 
>-* 0.5 
0.0 


G29.96-0.02 SH 



-20 0 20 40 60 80 100 

LSR velocity (km/s) 


G29.96-0.02 H 2 S 



G29.96-0.02 SH 



G29.96-0.02 HjS 



G29.96-0.02 CS 



G29.96-0.02 SO 




LSR velocity (km/s) 


G29.96-0.02 CS 




G29.96-0.02 HjS* 



LSR velocity (km/s) 


Fig. 5. Spectra obtained toward G29.96-0.02. Red sections of each plot indicate spectral regions assumed devoid of absorption for 
the purpose of fitting a continuum with emission lines. The right panels are zoomed versions of the left panels, with the scaling 
chosen to reveal weak absorption features. Magenta lines indicate the lambda doubling and hyperfine splittings for SH. 
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Table 4. Derived molecular column densities 


Source 

Vlsr 
(km s“') 

A(SH)‘> 

(10*2 cm-2) 

AfHzS) 

(10*2 cm-2) 

A(CS) 

(10*2 cm-2) 

A(SO) 

(10*2 cm-2) 

A(SH)/N(H2S) 

A(CS)/N(H2S) 

A(S0)/N(H2S) 

W49N 

37-44 

46.6(7.7) 

37.2(0.2) 

17.2(0.1) 

2.9 (0.3) 

1.25 

(0.21) 

0.46 (< 0.01) 

0.08 (< 0.01) 

W49N 

57-67 

22.1(8.3) 

13.0(0.2) 

6.9 (0.1) 

3.1 (0.4) 

1.70 

(0.63) 

0.53 (0.01) 

0.24 (0.03) 

W31C 

15-18 

34.5 (1.6) 

11.7(0.1) 

5.9 (0.1) 

2.5 (0.3) 

2.96 

(0.14) 

0.51 (0.01) 

0.21 (0.02) 

W31C 

18-21 

10.5(1.2) 

4.9 (0.1) 

1.6(0.1) 

0.6 (0.2) 

2.17 

(0.26) 

0.34 (0.02) 

0.12 (0.05) 

W31C 

21-25 

26.3(1.4) 

12.4 (0.2) 

4.9 (0.1) 

2.4 (0.3) 

2.13 

(0.12) 

0.39 (0.01) 

0.20 (0.02) 

W31C 

25-36 

113.9 (2.4) 

55.1(0.3) 

26.2(0.2) 

10.4(0.5) 

2.07 

(0.05) 

0.48 (< 0.01) 

0.19 (< 0.01) 

W31C 

36-45 

41.6(2.3) 

15.1(0.2) 

9.1 (0.2) 

4.5 (0.4) 

2.76 

(0.16) 

0.61 (0.02) 

0.30 (0.03) 

W51 

62 - 75* 

175.4 (2.1) 

155.5(1.7) 

37.0(0.5) 

63.3(1.1) 

1.13 

(0.02) 

0.24 (< 0.01) 

0.41 (< 0.01) 

G34.3-H0.1 

25-30 

5.9(1.1) 

4.3 (0.3) 

1.4(0.1) 

< 1.2(3 0-) 

1.36 

(0.28) 

0.32 (0.03) 

< 0.3 (3 cr) 

G29.96-0.02 

5-10 

21.8(4.1) 

19.5 (0.9) 

7.2 (0.2) 

1.6 (0.5) 

1.12 

(0.22) 

0.37 (0.02) 

0.08 (0.03) 


Numbers in parentheses are 1 sigma statistical errors 

® Material absorbing in this velocity interval may lie in close proximity to W51, where radiative excitation is important: the derived column 
densities, computed for the conditions in foreground diffuse clouds, should therefore be regarded with caution. 


Table 5. Derived molecular abundances 


Source 

Vlsr 
( km s“*) 

AlHz)” 

(10^* cm-2) 

A(H)* 
(10^* cm-2) 

MSH)/ 

A(H2)" 

X 10“* 

MHjS)/ 

A(H2) 

xlO"* 

N{CS)I 

N(U2) 

xlO** 

N{SO)/ 

N(H2) 

X 10** 

MSH)/ 

Wh 

xlO** 

MH 2 S)/ 

Wh 

xlO** 

MCS)/ 

A*h 

xlO"* 

N{SO)/ 

Nil 

X 10** 

W49N 

37-44 

1.76 

2.35 

26.5 

21.2 

9.8 

1.7 

7.94 

6.34 

2.92 

0.50 

W49N 

57-67 

4.19 

4.30 

5.3 

3.1 

1.7 

0.7 

1.74 

1.03 

0.55 

0.24 

W31C 

15-18 

2.34 

1.23 

14.7 

5.0 

2.5 

1.1 

5.84 

1.97 

1.00 

0.42 

W31C 

18-21 

1.57 

1.18 

6.7 

3.1 

1.0 

0.4 

2.44 

1.12 

0.38 

0.13 

W31C 

21-25 

2.49 

2.13 

10.6 

5.0 

2.0 

1.0 

3.71 

1.74 

0.68 

0.34 

W31C 

25-36 

8.19 

10.35 

13.9 

6.7 

3.2 

1.3 

4.26 

2.06 

0.98 

0.39 

W31C 

36-45 

4.94 

37.12 

8.4 

3.0 

1.9 

0.9 

0.89 

0.32 

0.19 

0.10 

G34.3+0.1 

25-30 

1.07 

3.72 

5.5 

4.0 

1.3 

< 1.1 {3a-) 

1.00 

0.73 

0.23 

<0.21 (3 cr) 

G29.96-0.02 

5-10 

1.43 

1.11 

15.3 

13.6 

5.1 

1.1 

5.50 

4.91 

1.82 

0.40 


Based on the CH column densities reported by Godard et al. (2012), for an assumed CH/H 2 ratio of 3.5 X 10“* (Sheffer et al. 2008), or - in the case of G29.96-0.02 
- on the HF column density reported by Sonnentrucker et al. (2015), for an assumed HF/H 2 ratio of 1.4 X 10“* 

® Obtained from 21 cm observations performed with the JVLA (Winkel et al. 2015) or the VLA (Fish et al. 2003, for the case of G29.96-0.02) 

The errors in the derived abundances are dominated by uncertainties in the CH/H 2 conversion factor (which is uncertain by a factor ~ 1.6) 


this assumption is justified by excitation calculations described 
in Appendix B. In the case of H 2 S, an ortho-to-para ratio of 3 
was assumed. Blue curves show multigaussian fits to the col¬ 
umn density per velocity interval. In the case of SH, the con¬ 
version to column density is complicated by the presence of 
hyperfine structure and lambda doubling. Here, the blue curve 
represents the fitted column density per unit velocity interval, 
while the red curve represents the corresponding optical depth 
(after convolution with the hyperfine structure and lambda dou¬ 
bling.) Thus, in this case, it is the red curve that should be com¬ 
pared with the optical depths shown in the black histogram. The 
blue symbols above the spectra indicate the centroids and widths 
(full-widths-at-half-maximum) for each fitted Gaussian. In our 
fits, the linewidths were constrained to be identical for all four 
molecules, and the centroids allowed to vary by at most 0.2 km/s. 

In Table 4, we present the foreground molecular column den¬ 
sities inferred for selected velocity intervals toward each source, 


along with column density ratios. The selected velocity windows 
are marked by the red vertical lines in Figures A.l - A.5, and 
the numbers appearing above each window give the molecular 
column densities in units of 10^^ cm“^. The average abundances 
for each velocity range appear in Table 5, computed relative to 
molecular hydrogen (columns 4-7) and relative to hydrogen 
nuclei (columns 8 - 11). For this purpose, we adopted atomic 
hydrogen column densities, A^(H), and hydrogen nucleus column 
densities, Ah = A(H) -H 2 A(H 2 ), obtained with the assumptions 
given in the table footnote. The values given in Table 4 for the 
W51 Vlsr = 62-75 km/s window must be regarded with cau¬ 
tion. Material absorbing in this velocity interval may lie in close 
proximity to W51, where radiative excitation is important: ex¬ 
cept for the case of SH, the derived column densities, computed 
for the conditions in foreground diffuse clouds, may be unreli¬ 
able. For this reason, the W51 Vlsr = 62-75 km/s has not been 
included in Table 5. 
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Fig. 6. Molecular column densities (top) and column density ratios (bottom) in foreground clouds along the sight-lines to W49N 
(blue), W31C (red), G34.3-1-0.1 (magenta) and G29.96-0.02 (black). The error bars indicate 1 cr statistical errors only. Dashed curves 
indicate the regions spanned by the range of column density ratios measured in this set of foreground clouds, with the dotted curve 
plotted midway between the dashed curves. The green point applies to the 62 - 75 km s * velocity range in W51, which probably 
lies in close proximity to the background continuum source. 


Our observations did not lead to any detection of H 3 S^. The 
line emission spectrum is quite rich in this spectral region, and 
the presence of interloper emission lines (i.e. lines from species 
other than H 3 S^) reduces our sensitivity to broad H 3 S^ fea¬ 
tures. In particular, an emission feature lying ~ 7kms * blue- 
ward of the H 3 S^ line - most likely the 32i - 4i4 transition 
of methanol at 293.464203 GHz - greatly limits our sensitiv¬ 
ity to H 3 S^ emission from the background sources. However, 
good upper limits can be placed on the H 3 S^ optical depths as¬ 
sociated with certain narrow H 2 S absorption components. The 
best limits on the optical depth ratio, t(H 3 S^)/t(H 2 S), are ob¬ 
tained for the 40 km/s absorption component in the foreground 
of W49N and the 30 km/s absorption component in the fore¬ 
ground of W31C. The 3 cr upper limits on t(H 3 S^)/t(H 2 S) ob¬ 
tained for these components are 1.06 x 10“^ and 1.02 x 10“^ 
respectively. Figure A.6 in Appendix A shows the H 3 S^ (green) 
and H 2 S (red, scaled by a factor 0.1) optical depth spectra ob¬ 
tained for these two components, together with the 3 cr upper 
limit on H 3 S^ (blue). These limits correspond to a column den¬ 
sity ratio A(H 3 S+)/A(H 2 S) ~ 0.028(7’,ot/100K). Here, is 
the rotational temperature for the metastable (7 = K) states of 


(the symmetric top molecule) H 3 S'^, a negligible population be¬ 
ing assumed for all other states within each /f-ladder. 

4. Discussion 

4. 1 . Measured column densities, abundances, and 
correlations among the molecules observed 

All the molecular abundances in Table 5 lie far below the solar 
abundance of sulphur, (ms/hh)© = 1.3 x 10“^ (Asplund et al. 
2009), indicating that the four molecules together account for 
at most 0.14% of the sulphur nuclei within any of the velocity 
ranges we defined. Moreover, sulphur is notable in showing no 
measurable depletion in diffuse atomic clouds (Sofia et al. 1994), 
and only moderate levels of sulphur depletion (by a factor 4) 
have been inferred even in the dense (uh ~ 10 ^ cm“^) gas within 
the Horsehead PDR (Goicoechea et al. 2006), suggesting that 
the molecules that we have observed contain a similarly small 
fraction of the gas-phase sulphur nuclei. 

Although the chemical pathways leading to the formation of 
SH, H 2 S, CS and SO are clearly inefficient in converting sul¬ 
phur atoms and ions into sulphur-bearing molecules, they evi¬ 
dently operate in a rather uniform manner, with the relative abun- 
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Fig. 7 . Molecular column densities (left) and abundances relative to H 2 (right) in foreground clouds along the sight-lines to W49N 
(blue), W31C (red), G34.3-1-0.1 (magenta) and G29.96-0.02 (black). The H 2 column densities are derived from CH column densities 
presented by Godard et al. (2012), for an assumed CH/H 2 ratio of 3.5 x 10“® (Sheffer et al. 2008). The error bars are dominated by 
uncertainties in the CH to H 2 conversion, except for the SO values toward G34.3-I-0.1. Dashed curves indicate the regions spanned 
by the range of abundances measured in this set of foreground clouds, with the dotted curve plotted midway between the dashed 
curves. 


dances among the sulphur-bearing molecules showing a remark¬ 
able constancy. In Figure 6, we show (top panels) the column 
densities of SH, CS and SO for each velocity interval in Table 3, 
as a function of N(H 2 S). The dependence is nearly linear in ev¬ 
ery case, with the column density ratios (bottom panels in Figure 
6 and columns 6 - 9 in Table 4) spanning a fairly narrow range; 


the A^(SH)/A^(H 2 S), A^(CS)/A^(H 2 S), and A^(SO)/A^(H 2 S) ratios 
lie in the ranges 1.12-2.96, 0.32 - 0.61, and 0.08 - 0.30 respec¬ 
tively. Furthermore, these column density ratios exhibit no obvi¬ 
ous trend with A^(H 2 S) as the latter varies over a factor ~ 13. In 
Figure 7, we present the column densities (left panel) and abun¬ 
dances (right panel) of SH, H 2 S, CS and SO relative to H 2 , as 
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Fig. 8. Molecular column densities (left) and abundances relative to H nuclei (right) in foreground clouds along the sight-lines to 
W49N (blue), W31C (red), G34.3H-0.1 (magenta) and G29.96-0.02 (black). The H 2 column densities are derived from CH column 
densities presented by Godard et al. (2012), for an assumed CH/H 2 ratio of 3.5 x 10“^ (Sheffer et al. 2008), and the atomic H column 
densities are derived from 21 cm observations reported by Winkel et al. (2015) or (for G29.96-0.02) by Fish et al. (2003). The error 
bars are dominated by uncertainties in the CH to H 2 conversion, except for the SO values toward G34.3-1-0.1. Dashed curves indicate 
the regions spanned by the range of abundances measured in this set of foreground clouds, with the dotted curve plotted midway 
between the dashed curves. 


a function of A^(H 2 ). Similar results are presented in Figure 8, 
but with abundances given relative to H nuclei as a function of 
A^h- Once again, there is no strong trend in abundance with either 

A^(H 2 ) or A^h- 


4 . 2 . Principal Component Analysis 

We have used Principal Component Analysis (PCA; e.g. Jolliffe 
2002 ) to explore differences and similarities among the absorp¬ 
tion spectra observed toward W31C and W49N. These sources 
are suited particularly well to this analysis because their absorp- 
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Fig. 9. Optical depths for eleven transitions observed toward 
W31C as a function of LSR velocity. Black curve: rebinned ob¬ 
served data (see text). Red histogram: approximate fit using only 
the first two principal components. Blue histogram: approximate 
fit using only the first three principal components. 

tion spectra exhibit a complex structure involving multiple fore¬ 
ground velocity components, detected at high signal-to-noise ra¬ 
tio, and unassociated with the source itself. In our analysis, we 
made use of the optical depths, as a function of LSR veloc¬ 
ity, for the SH, H 2 S, CS and SO transitions observed here, to¬ 
gether with analogous results for seven additional transitions of 
interest: the HI 21 cm transition (Winkel et al. 2015); the CH 
N - l,J - 3/2 - 1/2 lambda doublet (Gerin et al. (2010a); 
the H20^ 1 11 - Ooo transition (Gerin et al. 2010b; Neufeld et al. 


Fig. 10. Principal components obtained for W31C, from a joint 
analysis of the W31C and W49N absorption spectra 

2010); the ^^CH^ 7=1-0 transition (Godard et al. 2012); the 
SH^ Jn - I 2 - 0i transition; the ArH^ 7=1-0 transition 
(Schilke et al. 2014), and the CN N - I - 0 ttansition (Godard 
et al. 2010). In the case of those transitions (the SH, CH, SH^, 
'^CH^, and CN transitions) that show hyperfine structure, we 
used the hyperfine deconvolved spectra (i.e. the optical depths 
inferred for the strongest hyperfine component.) 

These eleven “optical depth spectra” were rebinned on a 
common velocity scale covering LSR velocities in the range 
10 to 50 km/s for W31C and 35 to 75 km/s for W49N. They 
are shown by the black histograms in Figure 9 (for W31C) and 
Figure C. 1 in Appendix C for W49N. Prior to implementing the 
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W31C and W49N 



W31C and W49N 



Fig. 11. Coefficients for the first two principal components ob¬ 
tained for each of the eleven transitions. Solid lines connect the 
origin to the point representing the coefficients. The black circles 
shows the unit circle within which all the coefficients must lie. 


PCA, we followed the standard procedure of scaling and shifting 
each optical depth spectrum so that the variances were unity and 
the means were zero. This ensures that each spectrum is given 
equal weight. The PCA method represents each optical depth 
spectrum as a linear combination of eleven orthogonal (i.e. un¬ 
correlated) components, with the optical depth for the /th transi¬ 
tion given by 

11 

Ti(\) ^ ai + biY^CijSj{v), ( 1 ) 

j=i 

where the S j(v) are the 11 principal component optical depth 
spectra shown in Figure 10 (for W31C) and Figure C.2 in 
Appendix C (for W49N), and a,, b,, and Cij are coefficients. The 


S j(y) all have a mean of zero and a variance of unity, and the Cij 
obey the normalization condition 2 / C?, = 1 ■ 

J IJ 

A key feature of PCA is that the terms within the sum in 
equation (1) above make successively smaller contributions to 
the right-hand-side, with the first term accounting for most of 
the variance in the optical depth spectra, the second term the 
next largest fraction of the variance, and so on. In this particu¬ 
lar case, the first two principal components are found to account 
for 72% of the total variance, and the red histograms in Figures 
9 and B.l show that a fairly good approximation to most of the 
observed spectra is obtained by retaining just the first two princi¬ 
pal component^ The only notable exception is ArH^, for which 
the first two terms in equation (1) yield a poor fit to the observed 
spectra; for that case, the blue histograms in Figures 9 and C.l 
show the sum of the first three terms. 

In Figure 11, we present, for each optical depth spectrum, /, 
the first three coefficients, C,i, Ca, and Ca in the sum in equa¬ 
tion (1). The upper panel shows Ca versus C,i, while the lower 
panel shows Ca versus Cq. Each transition occupies a specific 
position within these plots, which provides a valuable graphi¬ 
cal representation of the similarities and differences between the 
spectra. Color-coded squares show the coefficients for each tran¬ 
sition, and lines are drawn from the origin to each point. The 
results plotted in Figure 11 were obtained from a simultaneous 
analysis of the W31C and W49N data, but similar results were 
obtained by analysing each source individually. Moreover, the 
relative positions of the SH, H 2 S, CS, SO, CH, CN and HI tran¬ 
sitions in the upper panel do not depend strongly upon whether 
the molecular ion transitions are included in the analysis. Several 
features are noteworthy: 

First, with the notable exception of ArH^ transition, most 
transitions lie fairly close to the unit circle (black circle) in the 
upper panel, indicating that the first two principal components 
contain most of the information about their spectra. (If the coef¬ 
ficients for the third and higher principal components were zero, 
then all points would lie exactly on the unit circle.) ArH^ is un¬ 
usual in that its spectrum is largely accounted for by the third 
principal component. The fact that it is poorly correlated with 
most of the other molecules may reflect its unusual chemistry 
(Schilke et al. 2014): it is believed to be a tracer of gas with an 
extremely small ( i 10“^) molecular fraction. 

Second, the transitions of the four neutral sulphur-bearing 
molecules lie very close to each other in the plot, are signifi¬ 
cantly separated from the HI, CH, and all the molecular ion tran¬ 
sitions, and are relatively close to the CN transition. In the limit 
where the points lie very close to the unit circle, the correlation 
coefficient between any pair of transitions is equal to cosine of 
the angle between the relevant vectors (i.e. the lines drawn from 
the origin to each point). Indeed, the linear correlation coeffi¬ 
cients among the neutral sulphur-bearing molecules range from 
0.75 to 0.96 (with those among CS, H 2 S, and SH lying in the 
range 0.90 to 0.96). By contrast, the correlation coefficients be¬ 
tween CH and the neutral sulphur-bearing molecules are smaller, 
in the range 0.52 and 0.67, and those between HI and the neutral 
sulphur-bearing molecules are smaller yet: 0.32 to 0.39. 

Third, the molecular ions - other than ArH^ - have tran¬ 
sitions that typically lie close together and close to HI in the 
upper panel. In particular, the SH^ molecular ion shows a dis¬ 
tribution closer to that of CH^ and H 20 ^ than to that of the 
neutral sulphur-bearing molecules observed in the present study. 
The molecular ions do, however, occupy a range of positions in 

This feature of PCA, the ability to reduce the dimensionality of a 
dataset, is widely used in image compression. 
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Fig. 12 . Chemical network for sulphur- and carbon-bearing molecules (see text for details) 


the lower panel (i.e. they show a range of coefficients for the 
third principal component.) 

Fourth, the vector corresponding to CH in Figure 11 lies 
between the vectors for HI and those for CN and the neutral 
sulphur-bearing molecules. CH is believed to be a good tracer 
for H 2 , both in gas that is mainly atomic and in gas that is 
mainly molecular, while CN is thought to be abundant only in 
regions of large molecular fraction, so this behaviour may in¬ 
dicate that the neutral sulphur-bearing species are preferentially 
found in regions of large molecular fraction and are less preva¬ 
lent in partially-molecular regions where CH is present but CN 
is absent. We note, in this context, that the relative positions of 
CH^, CH, and CN in the upper panel of Figure 11 are consistent 
with previous UV absorption line observations (Lambert et al. 
1990; Pan et al. 2004, 2005), which showed that CH line profiles 
could be seen as the sum of a broad component similar to the 
CH^ profile, and a nari'ow one, similar to that of the CN line. 

4 . 3 . Comparison with previous observations of CS, H2S and 
SO 

Sulphur-bearing molecules have previously been studied within 
diffuse molecular clouds, both through observations of line emis¬ 
sion (e.g. Drdla et al. 1989; Heithausen et al. 1998; Turner 
1996a, 1996b) and absorption (Tieftrunk et al. 1994; Lucas & 
Liszt 2002; Destree et al. 2009). Absorption-line observations 
can provide the most reliable estimates of molecular column 
densities, because - as noted in Appendix B - the latter of¬ 
ten depend only weakly upon the assumed physical conditions 
in the absorbing material or upon the adopted rate coefficients 
for collisional excitation. Accordingly, we focus here on pre¬ 


vious absorption-line observations, which provide the most di¬ 
rect quantitative comparison with results obtained in the present 
study. 

Millimeter-wavelength observations of absorption by S- 
bearing molecules in diffuse clouds have been reported by Lucas 
& Liszt (2002), who detected CS, H 2 , and/or SO toward six 
extragalactic radio sources. The typical A(CS)/A(H 2 ) ratio in¬ 
ferred from these observations was 4x10“®, a value lying close to 
the middle of the range of A(CS)/A(H 2 ) obtained in the present 
study (Table 5). Absorption line observations of CS in diffuse 
clouds are also possible at ultraviolet wavelengths through ob¬ 
servations of the C-X (0,0) band near 1401A, which has been 
detected along 13 sight-lines toward hot stars (Destree et al. 
2009). The observed absorption line strengths are broadly con¬ 
sistent with the CS abundances inferred by Lucas & Liszt (2002) 
and in the present study, but a detailed comparison is limited by 
lack of knowledge of the relevant /-values. 

Lucas & Liszt (2002) reported N(CS)/N{SO) ratios in the 
range 1.7+0.8, which overlaps the range (2-6) implied by Table 
4, and A(CS)/A(H 2 S) ratios in the range 6+1. The H 2 S column 
densities derived by Lucas & Liszt (2002), however, appear to 
be based on a column-density-to-optical-depth ratio (their Table 
2 ) that is a factor of ~ 12 too small; applying this correction 
brings their A(CS)/A(H 2 S) ratio down to 0.48+0.08, in excellent 
agreement with the range (0.32 - 0.61) reported here in Table 4. 
In addition,Lucas & Liszt also detected the 85.348 GHz J - 2- 
1 transition of HCS^ in absorption toward one of the six sources, 
3C111, with an implied A(CS)/A(HCS+) ratio of 13.3+0.1. (The 
HCS^ J -2-1 frequency was not covered in the present study.) 
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Fig. 13. Column density ratios predicted by standard TDR mod¬ 
els, for a turbulent rate-of-strain of 3 x 10“" s“'; ion-neutral ve¬ 
locity drifts of 4 (blue), 6 (magenta), and 8 (red) kms^'; and 
densities of 30 (triangles), 50 (squares), 100 (circles), and 300 
(upside down triangles) cm“^. The column densities are com¬ 
puted for a total hydrogen column density A^h of 1.8 x 10 ^' cm“^, 
assuming an homogeneous shielding Ay over the entire line-of- 
sight. Along each curve. Ay varies between 0.1 and 1.0 from 
l|^t to right. Grey shaded regions indicate the observed range of 
values. 


4.4. Comparison with astrochemicai modeis 

4.4.1. Models without “warm chemistry” 

Standard models of photodissociation regions (PDRs) are 
known to greatly underpredict the abundances of sulfur-bearing 
molecules (e.g. Tieftrunk et al. 1994). For the densities (hh ~ 
few X 100 cm“^) and extinctions (Ay ~ few x 0.1 mag) typical 
of the diffuse foreground clouds we observed (e.g. Gerin et al. 
2014), we hnd the abundances to be underpredicted by a fac¬ 
tor of ~ 100, even in models that include gas-grain chemistry 
(Vasyunin & Herbst 2011) and the photodesorption and chemi¬ 
cal desorptiorO of molecules from grain surfaces. 

4.4.2. Models for turbulent dissipation regions 

As discussed by Neufeld et al. (2012), the discrepancy described 
in the previous subsection may provide evidence for a warm 
chemistry driven by turbulent dissipation in magnetized shocks 
(e.g. Pineau des Forets et al. 1986) or in vortices (TDR model, 
Godard et al. 2009) where the gas temperature and/or ion-neutral 
drift velocity are sufficient to activate highly endothermic for¬ 
mation pathways. Such a scenario has indeed proven success¬ 
ful (Godard et al. 2014) in reproducing the large abundances of 
CH^ and SH^ observed in the diffuse interstellar gas assuming 
that turbulence dissipates in regions with density tin < 50 cm“^, 
ion-neutral drift velocity Um ~ 3 km/s, and at an average rate 
e - 10“^"^ erg cm“^ s“', a value in agreement with measurements 
of the kinetic energy transfer rate at large scales (Hennebelle & 
Falgarone 2012). 

In Figure 12, we show the main production and destruction 
routes for 19 sulfur- and carbon-bearing molecules predicted 
by the PDR, TDR and shock models for a gas with density 
mh = 50 cm“^ and shielding Ay = 0.4 mag. For each species, 
only the processes that together contribute to more than 70% 
of the total destruction and formation rates are displayed. Solid 
blue arrows indicate reactions that are only important in PDRs 
while solid red arrows indicate reactions that are only impor¬ 
tant in TDRs or shocks. Green and dashed red arrows denote 
additional processes discussed in §4.3.4 below. Activation en¬ 
ergy barriers are indicated by temperatures (AEfk) marked next 
to the relevant reaction arrow. The chemical network we used 
here is an updated version of that used in the Meudon PDR 
code (Le Petit et al. 2006, http;//ism.obspm.fr). Table D.l in 
Appendix D lists the rates adopted for all the processes shown 
in Figure 12, together with those for two additional reactions for 
which updated rate coefficients have been used. Here, many of 
the rates adopted for endothermic reactions were obtained from 
those for the exothermic reverse reactions using the principle of 
detailed balance. This procedure may introduce inaccuracies, be¬ 
cause it is based upon the assumption that the reactants are in lo¬ 
cal thermodynamic equilibrium with the kinetic temperature. At 
the low densities typically of relevance in astrophysical applica¬ 
tions, the degree of vibrational excitation - in particular - may 
be highly subthermal. Thus, the reaction rates for endothermic 
reactions with vibrationally-cold reactants may be significantly 
overestimated if vibrational excitation is substantially more effi¬ 
cient than translational energy in overcoming the energy barrier. 
Unfortunately, detailed scattering calculations with a quantum 
chemistry potential are needed to determine which reactions are 
most affected; these are not typically available in the literature. 

^ Recent experiments by Dulieu et al. (2014) have suggested that the 
latter process is important in the case of water produced on bare grain 
surfaces. 
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TDR model predictions are shown in Figure 13, where the 
predicted column density ratios are plotted against A^(SH)/A^h- 
Grey shaded regions indicate the range of observed values listed 
in Tables 4 and 5 (or, in the case of A^(CS)/A^(HCS^), re¬ 
ported by Lucas & Liszt 2002). On the one hand, the observed 
A^(CS)/A^(SH) and A^(SH)/A^h ratios can be accounted for by 
models with an ion-neutral drift velocity of at least 6 km s“*, 
somewhat larger than that invoked to explain the CH^/SH^ col¬ 
umn density ratio. This result suggests that turbulent energy dis¬ 
sipates over a distribution of structures with different strengths 
rather than a single type of event. On the other hand, how¬ 
ever, the TDR model predictions systematically underestimate 
the column densities of both H 2 S and SO by roughly an or¬ 
der of magnitude unless Ay is unrealistically large (Ay > 1)0 
Moreover, the CS/HCS^ ratios reported by Lucas & Liszt (2002) 
for the sight-line toward 3C111 are larger than the model predic¬ 
tions except for the highest assumed density {tin = 300 cm~^). 

4.4.3. Shock models 

A similar discrepancy is found when the predictions of stan¬ 
dard shock models are compared with the data. Using the model 
described by Lesaffre et al. (2013), we have conducted a pa¬ 
rameter study in which the column densities of sulphur-bearing 
molecules behind a shock were computed as a function of the 
shock velocity, Vj, the preshock magnetic field, Bq, the preshock 
density, tin, and the UV shielding. The latter is parameterized 
by the visual extinction. Ay, to the cloud surface, which we 
assumed to be irradiated by the typical Galactic UV radiation 
field (Draine 1978). As long as Vj exceeds a value sufficient to 
drive the endothermic production reactions of relevance (about 
10 km/s for uh = 100 cm~^ and Ay = 0.1), the SH/H 2 S ra¬ 
tio is virtually independent of Vj and Bq. As in the TDR case, 
the predicted H 2 S and SO column densities relative to those of 
SH are underestimated by a factor of about 10 for clouds with 
hh = 100 cm“^ and small shielding Ay < 1. Similarly, the 
CS/HCS^ ratio is also underestimated except for high density 
values. Unlike TDR models, however, shock models underesti¬ 
mate the CS/SH ratio by a factor of about 5 for Ay < 1 . 

4.4.4. Additional effects on the chemistry 

As a possible explanation of the small observed V(SH)/V(H 2 S) 
ratios relative to the predictions of standard shock models, we 
have investigated two effects. First, we considered the depen¬ 
dence of the predicted SH/H 2 S column density ratio upon the 
branching ratio for dissociative recombination of HaS^. In the 
standard model, we adopted a branching ratio for H 2 S produc¬ 
tion of 0.17, the value measured by Kamiriska et al. (2008) for 
the production of D 2 S by dissociative recombination of D 3 S^. 
If we assume instead that every recombination of leads to 
the production of H 2 S, the SH/H 2 S ratio decreases to a value of 
about 5, i.e. about twice the typically-observed value. 

Second, we have investigated a physical effect not previously 
included in standard models for the chemistry of sulphur-bearing 
molecules in interstellar shock waves. After being produced by 
dissociative recombination of H 2 S^ or H 3 S^, SH molecules will 
initially carry an imprint of the velocity of the ionized parents 
from which they formed; this velocity drift, relative to other 


® This confirms the sense of the discrepancy discussed by Neufeld 
et al. (2012). However, Neufeld et al. (2012) overstated the magnitude 
of the discrepancy, because of a computational error that led to the ob¬ 
served SH/H 2 S ratio being understated by a factor 10. 



Fig. 14. Velocity of ions (red), SH (blue), and other neutral 
species (black), as a function of distance through the shock. 
Results are presented for two values of the rate coefficient, 
kmt(SH, H 2 ), for momentum transfer between SH and H 2 : 1 x 
10^^° (dashed blue curve) and 5 x 10“^°cm^s“' (solid blue 
curve). 


neutral species, might significantly enhance the production of 
H 2 S via the endothermic neutral-neutral reaction SH(H 2 , H)H 2 S 
(for which AEjk - 6862 K). As shown by Flower & Pineau 
des Forets (1998), who followed the velocity of CH molecules 
within a C-type shock, computing the drift between SH and H 2 
only requires the rate coefficient for momentum transfer between 
SH and H 2 , knit(SH, H 2 ). In principle, this rate could be obtained 
by subtracting the inelastic collisions rates (including both reac¬ 
tive and non reactive processes) from the total SH -H H 2 collision 
rate. However, as far as we know, the collision frequency for this 
system has never been estimated; moreover the rates of elastic 
and inelastic processes depend strongly on the gas kinetic tem¬ 
perature. As a result, we have explored the shock model predic¬ 
tions over a wide range of momentum transfer rate coefficient; 
10-'2 < kmt(SH, H 2 ) < 5 X 10- cm^ s-'. 

In Figure 14, we show the drift velocity of SH relative to H 2 , 
as a function of distance through the shock. As expected, the SH 
velocity lies somewhere between the ion velocity (red) and the 
H 2 velocity (black), with the exact behaviour depending upon 
the value adopted for kmt(SH, H 2 ). 

Figure 15 shows predictions for the shock model, with “ve¬ 
locity memory effects” included for SH. The reaction pathway 
enabled by such effects is shown by the green arrow in Figure 
12. The results are shown as a function of nn and for three values 
of Ay (0.1, 1, and 2). Here, the top panel shows the SH column 
density for a single shock, while the other four panels show the 
CS/SH, H 2 S/SH, SO/SH, and CS/HCS^ column density ratios, 
all computed for a single C-type shock with a velocity of 14 km/s 
and a preshock magnetic field Bo - (nH/cm-^)02^G. Results 
are shown for kmt(SH, H 2 ) of 1 X 10 '°cm^ s ' (dashed curves) 
and 5 x 10-'° cm^ s-* (solid curves). The grey shaded areas in 
each panel indicate the range of observed values presented in 
Table 4. In the case of the top panel, the shaded area refers to the 
SH column density per magnitude of visual extinction; a com¬ 
parison with the model predictions then leads to the number of 
shocks per magnitude of visual extinction required to reproduce 
the observations. For Ay < 1, the inclusion of velocity memory 
effects for SH, with an assumed knit(SH, H 2 ) of 1 x 10-*°cm^ s-* 
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Fig. 15. Column density predictions for the standard shock 
model, as a function of preshock density, hh, and for Ay of 0.1 
(blue), 1 (magenta), and 2 (red). Velocity memory effects are in¬ 
cluded for SH for assumed values of 1 x 10 '‘’cm^ s ' (dashed 
curves) and 5x10“'° cm^ s“' (solid curves) forkmt(SH, H 2 ). Gray 
shaded regions indicate the observed range of values. 


(dashed curves), significantly increases the predicted H 2 S/SH ra¬ 
tio, although the observed values still lie a factor ~ 3 above the 
predictions unless Ay is unrealistically large ( i 1). This process 
has however a small impact on all the other ratios displayed in 
Figure ISQ The SO/SH and CS/SH column density ratios are 
similarly underpredicted. The effect of velocity memory on the 
TDR model predictions will be investigated in future studies but 
is expected to have an impact similar to that found for shock 
models in increasing the predicted H 2 S/SH column density ra¬ 
tio. 

Finally, as a possible explanation of the underprediction of 
SO column densities by TDR and shock models, we have con¬ 
sidered a formation pathway of SO not yet included in stan¬ 
dard chemical networks and denoted by dashed red arrows in 
Figure 12: H 2 S^ -H O ^ HSO^ + H followed by the disso¬ 
ciative recombination of HSO'^. So far, standard chemical net¬ 
works available online (e.g. the UMIST database; McElroy et 
al. 2013) only include the channels H 2 S^ -H O ^ SH^ + OH 
and H 2 S^ -h O ^ SO^ + H 2 with identical rate coefficients of 
3.1 X 10“'°cm^s“' (Prasad & Hunttess 1980), despite the fact 
that the production of HSO^ is also a highly exothermic process. 
Implementing this channel in the shock model with an assumed 
rate coefficient of 6 . 2 x 10 “'° cm^ s ' enhances the predicted col¬ 
umn density of SO by a factor of 10. Such a scenario has yet to 
be confirmed by theoretical or experimental estimates of the re¬ 
activity of H 2 S^ with oxygen atoms and a detailed analysis of 
the possible products. 

4.4.5. Line widths and velocity shifts 

In the case of TDRs, an interpretation of the observed line pro¬ 
files requires a convolution of the kinematic signatures of large 
scale turbulence with those of active regions of dissipation - 
taking into account their orientations, strengths and distribution 
along the line of sight - and of the relaxation stages that nec¬ 
essarily follow the dissipative bursts. In TDRs, SH^, CH^, and 
H 3 S+ are formed in the active phase (Godard et al. 2014), while 
SH, H 2 S, SO, and CS are mainly produced by the relaxation 
phase. The linewidths of the former species are therefore ex¬ 
pected to be wider than those of the latter; this effect is indeed 
apparent in the optical depth spectra (Figures 9 and C.l) and 
shows up in the PCA results (Figure 11). In C-type shocks, the 
chemical production may induce a velocity shift between the line 
profiles of species. Flower & Pineau des Forets (1998) estimate 
the shift between CH and CH^ induced by a single shock of 14 
km/s propagating along the line of sight to be roughly 2 km/s. 
The same could be expected between SH and the neutral fluid, 
yet no significant line shifts are apparent in the observed data. 
We note, however, that taking a random distribution of shocks 
with different orientations will drastically reduce the line shifts. 
A detailed calculation of the line profiles in the framework of 
both the TDR and shock models remains to be done. 


5. Summary 

1. With the use of the GREAT instrument on SOEIA, we 
observed the SH 1383 GHz ^Uy 2 J = 5/2 - 3/2 lambda 
doublet towards the star-forming regions W31C, G29.96-0.02, 


’’ Predictions obtained for very small values of k[„t(SH, H 2 ) smaller 
than 3 x 10“" cm^ s“' (not shown) are in agreement with the observed 
H 2 S/SH column density ratio even for very diffuse {iin < 50cm“^) and 
weakly shielded (Ay = 0.1) environments. 
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G34.3-I-0.1, W49N and W51. We detected foreground absorp¬ 
tion toward all five sources. 

2. Using the IRAM 30m telescope at Pico Veleta, we observed 
the H 2 S 1 10 - loi (169 GHz), CS 7 = 2 - 1 (98 GHz), SO 32 - 2i 
(99 GHz) and HsS^ Iq - Oq (293 GHz) transitions toward the 
same five sources. The H 3 S^ Iq - Oq transition was not detected 
in any source, but, for the other three transitions, we detected 
emission from the background source, and absorption by fore¬ 
ground clouds, toward every source. 

3. In nine foreground absorption components detected to¬ 
wards these sources, the inferred column densities of the 
four detected molecules showed relatively constant ratios, with 
N(SH)/N(H 2 S) in the range 1.1 - 3.0, A(CS)/A(H 2 S) in the 
range 0.32 - 0.61, and A(SO)/A(H 2 S) in the range 0.08 - 0.30. 
The observed SH/H 2 ratios - in the range 5 to 26 xlO ® - indi¬ 
cate that SH (and other sulphur-bearing molecules) account for 
«; 1% of the gas-phase sulphur nuclei. 

4. We have used Principal Component Analysis to investigate 
similarities and differences in the distributions of SH, H 2 S, CS, 
SO, and seven other species present in foreground diffuse clouds 
along the sight-lines to W31C and W49N. The neutral sulphur¬ 
bearing molecules we observed are very well correlated amongst 
themselves and with CN, moderately well correlated with CH 
(a surrogate tracer for H 2 ), and poorly correlated with atomic 
hydrogen and the molecular ions SH^, CH^, H20^, and ArH^. A 
natural, but not unique, interpretation of these correlations is that 
SH, H 2 S, CS, SO are preferentially found in material of large 
molecular fraction. 

5. The observed abundances of neutral sulphur-bearing 
molecules greatly exceed those predicted by standard models of 
cold diffuse molecular clouds, providing further evidence for the 
enhancement of endothermic reaction rates by elevated temper¬ 
atures or ion-neutral drift. 

6 . Standard models for C-type shocks or turbulent dissipation 
regions (TDRs) can successfully account for the observed abun¬ 
dances of SH and CS in diffuse clouds, but they fail to reproduce 
the column densities of H 2 S and SO by a factor of about 10. 
More elaborate shock models, which include “velocity mem¬ 
ory effects” for SH, reduce this discrepancy to a factor ~ 3. 
Moreover, a possible reaction pathway not included in previous 
chemical networks, H2S^(0, H)HSO^(e, H)SO, can potentially 
explain the larger-than-expected SO abundances. 

7. Further theoretical studies will be needed to determine the im¬ 
portance of velocity memory effects in TDRs and any kinematic 
signatures that might be expected. 
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Appendix A: Absorption features and moiecuiar 
coiumn densities 

Figures A.l through A.5 show the column density per veloc¬ 
ity interval for each detected molecule in each observed source. 
Figure A.6 presents upper limits on HaS^. 

Appendix B: Excitation of suiphur-bearing 
moiecuies in diffuse moiecuiar ciouds 

In order to determine the equilibrium level populations for CS 
and SO in diffuse molecular clouds, we have solved the equa¬ 
tions of statistical equilibrium as a function of temperature and 
density. Our excitation model includes radiative pumping by the 
cosmic microwave background (CMB), at temperature Fcmb = 
2.73 K, spontaneous radiative decay, and collisional excitation 
by molecular hydrogen. We adopted the collisional rate coef¬ 
ficients recommended in the LAMDA database (Schoier et al. 
2005j3, which were obtained by a simple scaling of the values 
computed for the excitation of CS (Lique et al. 2006a) and SO 
(Lique et al. 2006b) by He. 

In Figure B.l, we present results for the ratio of total molec¬ 
ular column density to optical depth, NIt, for the CS 7 = 2 - 1 
(upper panel) and SO Jn -'ii- 2i (lower panel) transitions that 
we observed. Results are presented as a function of H 2 density 
for four different temperatures. The left panels show the results 
on a log-log plot, while the right panels show the same data on 
a log-linear plot for a narrower range of n(H 2 ). In the limit of 
low density, the level populations are in equilibrium with the 
CMB, at a temperature of 2.73 K, while in the limit of high 
density the level populations are in local thermodynamic equi¬ 
librium (LTE) at the gas temperature, T. In the latter limit, A/t 
is proportional to provided kT » the transition energy, AF ; 
this result arises because the difference between the absorption 
and stimulated emission rates is proportional to 

(Ni-Nugilgu) ^ A,[l-exp(-A£/kr)] ex A,/r (x A/(Z7’) (x T-\ 

where Nu and A; are the column densities in the lower and up¬ 
per states, and Z is the partition function, proportional to T for 
diatomic molecules. At intermediate densities and temperatures 
i 100 K, the populations are inverted. 

Two caveats affecting the results shown in B.l are our ne¬ 
glect of radiative trapping and of electron impact excitation in 
determining the level populations. While the SO transition is 
clearly optically-thin in all the absorption components we ob¬ 
served, the optical depth of the CS transition can reach values 
of order unity at line center; thus, for CS, radiative trapping can 
lead to a modest reduction in the critical density at which de¬ 
partures from the low-density behaviour become significant. For 
polar molecules like CS and SO, electron impact excitation can 
be of comparable importance to excitation by neutral species in 
regions where carbon is fully-ionized. If the CS and SO absorp¬ 
tion is occurring in regions where C is ionized, a further modest 
reduction in the critical density could result. 

Notwithstanding these caveats, for the temperatures {T < 
80 K) and densities (< few x 10^ cm“^) typical of the foreground 
diffuse molecular clouds observed in this study (e.g. Gerin et 
al. 2014), the A/t ratio is very comfortably in the low density 
limit for both CS and SO, justifying the assumption underlying 
our computation of column densities in §3 above. In the case 
of H 2 S, the effects of collisional excitation upon the A/t ratio 


http://home.strw.leidenuniv.nl/ ~ moldata/ 


are expected to be even smaller than those for CS and SO, both 
because the spontaneous radiative rate for the H 2 S transition is 
larger than those for CS and SO, and because the H 2 S absorp¬ 
tion takes place out of the lowest rotational state of the ortho 
spin-species. 

Appendix C: Principal component analysis for 
W49N 

Figures C.l and C.2 show the optical depth spectra for eleven 
transitions observed towards W49N, and the principal compo¬ 
nents determined from a simultaneous analysis of the W31C and 
W49N spectra. 

Appendix D: Reaction network for sulphur-bearing 
molecules 

The chemical network for sulphur-containing molecules adopted 
in the Meudon PDR code version 1.4.4 - a slightly-updated ver¬ 
sion of which is used here for TDR and shock modeling - in¬ 
cludes 432 reactions. In Table D.l, we present an abbreviated 
list containing rates for those reactions and protoprocesses that 
play a dominant role in the formation and destruction of sulfur 
bearing species (Figure 12) and/or whose rates have been up¬ 
dated for use in our study. 
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W49N 






LSR velocity (km/s) 


Fig.A.l. Molecular column densities along the sight-line to W49N. Black histograms show the data and blue curves show a multi- 
Gaussian fit to the data. In the case of the SH data, the red curve indicates a convolution of the multi-Gaussian fit with the hyperfine 
structure; it is this curve that should be compared with the black histogram (see text). Red boxes indicate the velocity ranges for 
which results are presented in Tables 1 and 2, with the red numbers at the top showing the integrated column density in units of 
10'^ cm“^. Blue symbols above each panel indicate the velocity centroid and FWHM for each Gaussian absorption component. 
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W31C 






Fig. A.2. Same as Figure A.l, but for W31C. 
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W51 



LSR velocity (km/s) 
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Fig.A.3. Same as Figure A.l, but for W51. 
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Fig. A.4. Same as Figure A.l, but for G34.3-1-0.1. 
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G29.96-0.02 






Fig. A.5. Same as Figure A.l, but for G29.96-0.02. 
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W49N 



LSR velocity / km s 


W31C 



LSR velocity / km s 


Fig.A.6. H 3 S^ (green) and H 2 S (red, scaled by a factor 0.1) optical depth spectra obtained for the W49N 40 km/s foreground 
absorption component and the 30 km/s W31C foreground absorption component. The 3 cr upper limit on is shown in blue. 
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Iog,o N(SO) in cm ^ per km/s for t=1 log,o N(CS) in cm ^ per km/s for 
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CS J=2-1 



SO Jfg = 32-2, 



CS J = 2-1 



log,o n(H2)/cm 



Fig. B.l. Column densities required for an optical depth of unity. 
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W49N SH 




W49N CS 



W49N SO 



S, for W49N 





S 4 for W49N 



W49N HI 



W49N CH 





W49N CN 



W49N 






Fig. C.l. Optical depths for eleven transitions observed toward Fig. C.2. Principal components obtained for W49N, from a joint 

W49N as a function of LSR velocity. Black curve; rebinned ob- analysis of the W31C and W49N absorption spectra 

served data (see text). Red histogram; approximate fit using only 

the first two principal components. Blue histogram; approximate 

fit using only the first three principal components 
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Table D.l. Abbreviated reaction list 


Reaction" 

A 

B 

c 

Reference" 

S + Hz ^ SH -t H 

4.20 X 10^'° 

0.00 

10494.00 

( 1 )“' 

SH + Hz ^ HzS -H H 

7.36 X 10-‘3 

2.31 

6862.00 

( 2 ) 

SH -t- H -> S Hz 

1.39 X 10-'° 

-0.18 

697.00 

( 1 ) 

S + OH ^ SO + H 

6.59 X 10-" 

0.00 

0.00 

(3) 

SH -t 0 SO -t H 

1.58 X 10-‘° 

0.00 

0.00 

(3;4) 

S+ -I- Hz ^ SH+ -H H 

2.20 X 10-‘° 

0.00 

9860.00 

(5)“' 

S+ -H Hz ^ HzS+ -H hv 

1.00 X 10-'^ 

- 0.20 

0.00 

( 6 ) 

SH+ -(- Hz ^ HzS+ -H H 

1.90X 10-‘° 

0.00 

6380.00 

(5)“' 

SH+ -l-Hz ^ H3S+ -Hhv 

2.40 X 10-'® 

-0.80 

0.00 

(7) 

SH+ H S+ -H Hz 

1.10 X 10-'° 

0.00 

0.00 

(5) 

HzS+ -H Hz ^ H 3 S+ + H 

1.40 X 10-" 

0.00 

2900.00 

(5)" 

HzS+ -r H ^ SH+ -(- Hz 

2.00 X 10-‘° 

0.00 

0.00 

(5) 

H3S+ -hH ^ HzS+ -hHz 

6.00 X 10-" 

0.00 

0.00 

(5) 

CS+ -H Hz ^ HCS+ - 1 - H 

4.50X 10-‘° 

0.00 

0.00 

( 8 ) 

S+ -H CH ^ CS+ -H H 

6.20 X 10-'° 

- 0.00 

0.00 

{9Y 

SH+ -H C ^ CS+ + H 

9.90 X 10-'° 

0.00 

0.00 

(9) 

CH+ -t- S ^ HCS+ -H Hz 

1.40 X 10-° 

0.00 

0.00 

(9) 

H+ -H S ^ SH+ -H Hz 

2.60 X 10-° 

0.00 

0.00 

(9) 

SH+ -H e ^ S -H H 

2.00 X 10-’ 

-0.50 

0.00 

(9) 

HzS+ -H e ^ SH -t H 

2.40 X 10-’ 

-0.72 

0.00 

( 10 )"' 

HzS+ -te^S-t-H-t-H 

2.40 X 10-’ 

-0.72 

0.00 

( 10 )"' 

HzS+ -H e ^ HzS -H hv 

1.10 X 10-'° 

-0.70 

0.00 

(9) 

H 3 S+ -re ^ SH-hH + H 

3.02 X 10-’ 

-0.50 

0.00 

( 11 ) 

H3S+ -re ^ HzS -t-H 

8.84 X 10-" 

-0.50 

0.00 

( 11 ; 12 ) 

H3S+ -re ^ SH-t-Hz 

7.80 X 10-" 

-0.50 

0.00 

( 11 ; 12 ) 

H 3 S+ -re ^ S -i-Hz -tH 

5.20 X 10-" 

-0.50 

0.00 

( 11 ; 12 ) 

HCS+ -)- e ^ CH -H S 

7.86 X 10-’ 

-0.57 

0.00 

(13) 

HCS+ -t e ^ CS -1- H 

1.84 X 10-’ 

-0.57 

0.00 

(13) 

0 -t- HzS ^ SH + OH 

7.90X 10-‘2 

1.00 

1737.00 

( 1 )"' 

H+ -t-HzS ^ H 3 S+ -hHz 

3.70 X 10-° 

0.00 

0.00 

(14;15)‘' 

Photoprocess*” 

(o 

r 


Reference" 


S + hv S+ -t e 

4.40 X 10-'° 

2.58 

(16) 

SH -t- hv ^ S -t- H 

6.50 X 10-'° 

1.42 

( 16 ) 

HzS -H hv ^ SH + H 

1.00 X 10 ° 

1.87 

(16;17) 

HzS -H hv ^ S -H Hz 

1.00 X 10-° 

1.87 

(16;17) 

CS - 1 - hv ^ C -H S 

6.30 X 10-'° 

2.03 

( 16 ) 

SO -t hv ^ S -H 0 

2.40 X 10-° 

1.95 

( 16 ) 

HCS+ -t hv ^ CS+ -H H 

7.00 X 10-'° 

1.70 

(18) 


Reaction rate coefficient k{T) - A(T /300K)® exp(-C K/T) cm^ 

® Rate if = (To expC-yAv/mag) 

References: (1) Manion et al. (2013); (2) Peng et al. (1999); (3) Demore et al. (1997); (4) Tsuchiya et al. (1994); (5) Millar et al. (1986); 
( 6 ) Herbst et al. (1989); (7) Millar & Herbst (1990); ( 8 ) Decker et al. (2001); (9) Prasad & Huntress (1980); (10) Hellberg et al. 2005; 
(11) Abouelaziz et al. (1992); (12) Kaminska et al. (2008); (13) Montaigne et al. (2005); (14) Kim et al. (1974); (15) Anicich et al. (1975); 
(16) van Dishoeck et al. (1988); (17) van Dishoeck (2006); (18) estimate 

Value different from that adopted previously in the Meudon PDR code version 1.4.4 
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